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numericai modelling free of experimental uncertainties. It was then applied to the Naval 
Postgraduade School axial research compressor with which tip clearance effects were to 
be investigated experimentally. The evaluation identified both limitations in the code 
and hmitations in the data which could be obtained in the experiment. It was 
recommended that provision be made to accommodate peripherally non-uniform flow 
effects (blockage) in the throughflow code calculation, and that provision be made to 


obtain adequate peripheral flow survevs in the experiment. 
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I. INTRODUCTION 


The effect of finite rotor tip clearance on the operating characteristics of 
compressors is not well understood. An investigation is currently underway at the 
Turbopropulsion Laboratory of the Naval Postgraduate School (NPS) to explore the 
interaction of the tip leakage flow with the rotor passage flow, and to identify 
parameters which control the effect of tip clearance on compressor performance. 
Experiments are conducted on a 36 inch diameter low speed axial flow compressor test 
rig. The compressor has three repeating stages and is currently configured with two 
stages of symmetrical blading. Measurements of the compressor performance are 
obtained using fixed probes and rakes, a flow nozzle and a torque meter. Flow survey 
data are obtained beween blade rows using pneumiatic probes. Rotor exit conditions are 
measured using hot-wire probes. Also, using high-response semi-conductor pressure 
transducers and conditional sampling, the pressure field on the axisymmetric surface 
across the rotor tip is mapped. A complete description of the three dimensional flowfield 
within the compressor is desired. H1owever, for each type of flow survey, the coverage 
of the compressor annulus is limited by the access holes provided in the heavy 
compressor case-wall. Since the experimental program would be expected to develop a 
modification to current models for tip-clearance effects on throughflow development, a 
computational code incorporating a case-wall boundary layer calculation was needed 
with which to obtain predictions of the experimental conditions. The code would then 
provide a vehicle in which to incorporate a change in the tip-flow model. 

The finite element code Q3DFLO’81 was examined with this purpose in mind. The 
code was installed on the NPS mainframe IBM 370-3033 computer in 1983 by its 
Originator, Professor Charles Hirsch, and was used sucessfully at that time. 

The purpose of the present evaluation was specifically to determine the suitability 
of the code for the tip-clearance application, for which the case-wall boundary laver 
modelling would be a significant consideration. The evaluation was performed in two 
steps. First, the code was applied to a fully documented, ‘benchmark’, test case wherein 
completely detailed flow field information was available in addition to the geometry and 
the controlled boundary conditions. The findings from this effort are described in 
Chapter III. Then, with a developed appreciation of the degree to which the code 


described actual compressor flow conditions, the code was applied to the NPS 


compressor geometry. A comparison was made of code predictions with the 
pressure-rise performance measured experimentally. The results are described in 
Chapter IV. In order to provide a comparison with measured rotor-tip pressure 
contours, blade-to-blade calculations were also made at two compressor throughflow 
rates. An example 1s described in Chapter V. Conclusions and recommendations 
following both phases of the evaluation are given in Chapter VI. In general, it was 
found that the most significant limitation of the code for the desired application hes in 
the absence of a mechanism to input peripherally non-uniform flow effects, which are 
certainly always present in the machine. The need to obtain complete peripheral flow 
surveys in the experimental program was also identified. 

Before detailing the results of the evaluation, a review of the theoretical background 


and modeling on which the code 1s based, 1s given in Chapter I. 


Il. PRINCIPLES OF TURBOMACHINERY THROUGHFLOW 
MODELLING AND DESCRIPTION OF THE CODE 


A. BACKGROUND 

Q3DFLO’81 is a Finite Element code developed by Prof Charles Hirsch and his 
coworkers at Vrye University in Brussels for computing the complex flowfield in a 
turbomachine and its associated axisymmetric ducting. The formulation of the 
computational model of the flowfield to its final form for numerical coding is described 
in references | and 2. Reference | describes all available approaches to the treatment 
of turbomachinerv flowfields and derivative computational methods, including those 
Meedn tlie development of Q3DFLO’81. Reference 2 1s a summary of the key equations 
programmed in the Q3DFLO’81 computational code, which was the focus of the present 
study. 

Since the two references were intended for specialists in the field of turbomachinerv 
flow modelling, no detailed algebraic operations Were provided from one step to another. 
For a newcomer like the author, it was not easy to follow the logic and the mathematical 
language used in the presentation. To complicate matters, the equation of motion can 
be written in several different forms, although all are equivalent. Tvpicallv. for a flow 
in an annulus region, the equation is normally expressed with respect to an absolute 
frame of reference and within the blade to blade region with respect to a relative frame 
of reference. However, in the Q3DFLO‘’81 formulation, the equation of motion contains 
variables from both the absolute frame and the relative frame, in the same equation! 

The objective of this chapter therefore is to review the development of the equations 
which are the basis of flowfield computations in turbomachinery. The relationship 
between the various forms of the equation of mouon is shown and the detailed algebraic 
Operations that were omitted in summarizing the development of the Q3DFLO‘’81 in 
reference 2 are described here. Detailed derivation of equations which are available in 
texts and references will not be repeated. In such situations, only the key equations and 
their physical interpretations are presented. It is hoped that this will give simular 
newcomers a basic understanding before attempting the advanced materials in 


reference |. 


B. NON-AXISYMMETRIC NATURE OF THE REAL FLOW 

Figure | shows a blade row with radial leading and trailing edges and 4 families of 
surfaces; namely, surfaces of revolution S, , blade surfaces S, and stream surfaces S, and 
S,. S, is the axisymmetric surface generated by a generatrix m, as it rotates 360° about 
the axis. S, 1s simply the blade surface. Suppose a row of dyed particles is initially held 
along the edge A-P-C and another row along B,-P-B, at the inlet. When these particles 


are released, 1b 1S foumel tna 


e The row of particles initially along A-P-C does not follow the blade surface Sz. 
They instead follow the stream surface S, contrary to one’s expectation. 


e Similarly, the row of particles initially along B,-P-B, does not follow the 
axisymmetric surface S,. Instead they follow the stream surface S,. 
These observations have been verified by experiments. A flow through stationary or 
rotating cascades cannot be axisymmetric if the flow is to exert a moment on the blades. 
PRR eek 2186) 

In the annulus region between the blade rows, the flow can be considered to be 
axisymmetric after the wake has mixed out. In the formulation of through flow analvsis, 
the flow within the blade row is also very often assumed to be axisymmetric to simplify 
the mathematical model 1.e. surfaces S. will coincide with S, and S, Withee 
assumption 1s also adopted in the following derivation. However, this assumption 
causes discrepancies in the velocity components and flow angles between experimental 
measurement and throughflow modelling results. These discrepancies can be corrected 
by incorporating a factor known as the tangential blockage. This is dealt with in the 


next chapter. 


C. EQUATION OF MOTION IN THE ABSOLUTE FRAME 


The fundamental equation of fluid motion can be expressed as 


Ay VP 





A i = 
at Ve VV =——5-- WZ) + Fy {1} 


for viscous, compressible or incompressible fluids. 
In turbomachinery applications, the equation of motion is very often expressed in 
terms of enthalpy instead of pressure. This is done in 2 steps. Furst, use the vector 


identity 


VN SVS eae {2} 


Second, introduce the relationship between two thermodynamic states along the path 


of a fluid particle 








Vh— TVS =76 (3} 
Using Eq. (2) and Eq. (3), Eq. (1) becomes 
eV a = 2 
ee ome ae IVS Vi Vig Z) + Fy {4} 
Re-arranging, 
aT? 72 a a 
J+ V(h+ H+ gZ)=TVS+VKVXV +E, (5} 
The total enthalpy is defined as 
Vy? 
H=h+ 5 es {6} 
Reine Eq. (6), Eq. (5) becomes 
AL e een 
+ VH=TVS+VKVKV +E, (7} 


iiieedetailed derivation of Eqn (2) and (3) is given in reference 3. For a steady and 


isentropic flow, Eq. (7) reduces to 
VH=VxVxV {8} 


If V7 = 0, the fluid is either non-rotational or the vorticity is parallel to the particle 


velocity. 


D. EQUATION OF MOTION WITHIN A BLADE ROW 

Figure 2 shows the outlet flow of a constant speed rotating rotor with a relative 
velocity profile represented by the vectors HV, to IV,. A stationary observer attached to 
the stator blade row looking at the rotor outlet will see a non-steady velocity profile 
represented by the vectors V, to V, which vary in both magnitude and direction. It is 


clear that Eq. (7) is not useful for describing the flow within the rotor region, as the term 


av 





cannot be neglected nor determined. Hence, it 1s more useful to adopt a rotating 
frame of reference that rotates at the same speed as the rotor for describing the flow 


Within the rotor blade row. [Ref. 3: pp. 110,111] 


E. EQUATION OF MOTION IN THE RELATIVE FRAME 
The equation of motion with respect to the rotating frame can be obtained directly 


from Eq. (1) by replacing the absolute velocity by the relative velocity using 


—> —* 


V = aor {9} 





in place of re The left hand side (LHS) of Eqn (1) can be expressed as 














ay; = am = 

= +7 Vi =<(W +x R)+ (0 + Ox R) VV +x R) 
{10} 

op ~ ae 
= ~ aii? Vell’ +20 x W— V al ae ) 
and Eq. (1) becomes 
AV 2 °R? Mee 
a+ -v pil 20m i Vig ae ) = - EE _ va(gZ) + F {11} 


The detailed derivation of Eqn (10) 1s given in reference 3 . 
Again using the vector identity of Eqn (2) and the thermodynamic relationship of 


Eqn (3) in the relative frame, Eqn (11) becomes 
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{12} 
= TVpS — Vah—ValeZ) + F 


Re-arranging, 


al 2 R? 


Cal 








)=TVpS+Wx(VpxW+20)+F, {13} 


With the rothalpy defined as 


yo 2 pn2 
[=h+— + gZ- we {14} 





the equation of motion in the relative frame is then 


a, 
Ct 





+ Val = TVp5 + Wx (Vp x W + 20) + F, {15} 


For a steady and isentropic flow within the rotor and with the rothalpy constant, Eq. 
Geo ereduces to 


IV x (Vax W+20)=0 {16} 


Physicallv, this means that the flow is rotational with respect to the rotating frame and 
the relative vorticity 1s opposite to and equal to twice the rotor rotation. 

However, it should be pointed out that the assumption of steadv relative flow, 
though frequently made in turbomachinery calculations, is true only if the inlet flow to 
the rotor is uniform in the tangential direction: 1.e. the wakes of anv upstream blade rows 
must have mixed out. Hirsch [Ref. 1: pp. 53] states that in a real flow, “the viscous 
interactions and unsteadiness, which are always present to some degree in 


turbomachinery flow, introduce variations in rothalpy along the flow path.” 


F. EQUATION OF MOTION USED IN Q3DFLO’81 
An alternative form of the equation of motion for steady, relative, three-dimensional 
flow which is used in Q3DFLO’81 involves both relative and absolute velocity 


components, namely, 


—> ——b 


Val = 1VpS + Wx (Vx V)4 FF a7} 


foamy) 1s the same as Eq. (15) for a steady flow, except that the term (V, x Hr + 2w) 


is expressed in terms of (V x V). 


G. SIMPLIFIED THROUGHFLOW MODEL OF Q3DFLO’81 

Figure 3 shows the characteristics of an actual flow within the blade passage of a 
rotor. The flow is highly three-dimensional with strong corner vortices generated at the 
corners between the blade and hub and between the blade and shroud. There is also a 
secondary flow surrounding the primary core flow. However, the Q3DFLO’81 


throughflow uses a simplified axisymmetric flow model, and the three-dimensional 


details are removed (accounted for) by taking the density average of the flow over the 
blade spacing at each spanwise position. 

In assuming the flow to be axisvmmetric through blade row, the effects of the blades 
themselves must be introduced into the equation of motion as a body force. The 


equation of motion (Eq. 17) becomes 


ome 


—Wx(Vx V) = 1V,5 — Vere F, + F; {18} 


H. BASIC FORMULATION OF Q3DFLO’81 

The basic steps of the formulation of Q3DFLO’81 throughflow model are described 
in reference 2. The throughflow equation known as the Radial Equilibrium Equation 
that is solved in the Q3DFLO’81 code for the flow within the rotor is 
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(19) 1s derived from Eq. (18). Details of the derivation are given in Appendix (A). 

For the flow within the stator, Eq. (19) is still applicable. As the stator 1s stationary, 
the relative velocity components are replaced by the absolute velocity components and 
the rothalpv 1s replaced by the enthalpv. 

However, for the flow within the annular region without blades such as the annular 
space between a rotor and stator blade rows, Eq. (19) is not applicable. The inviscid 
equation of motion is used in such regions. The formulation of the form of equation 
of motion that is solved in Q3DFLO’81 for the flow in blade-free annular regions 1s 
given in reference 2. 

Figure 4 shows the boundary conditions applied to the throughflow calculation. 
The last station must be be normal to the meriodional velocity. 

The radial equilibrium equations for the flow within the blade row and the annular 
region are non-linear. They are solved by the finite element method using a iterative 


procedure. The details of the numerical method are also described in reference 2. 


II. ASSESSMENT OF THE CODE USING A BENCHMARK DATA SET 


A. BACKGROUND 

The set of radial equilibrium equations, which is solved in Q3DFLO’S1 is only one 
of three major ingredients in the code. Q3DFLO’81 incorporates an endwall boundary 
layer calculation which (optionally) can be coupled into the radial equilibrium solution 
procedure. The flow is then solved iteratively to account for the effect of boundary layer 
and tip clearance losses. Like all other throughflow codes, Q3DFLO’81 also makes use 
of correlations of two dimensional cascade data for the loss coefficient and flow turning 
angle. Consequently, the ability to predict the performance of an axial compressor using 
Q3DFLO’81 depends not only on the simplified modelling of the physical flow and the 
accuracy of the numerical method but also on the applicability and accuracy of the two 
dimensional cascade data. 

Comparisons of predictions using Q3DFLO’81 with experimental data on axial flow 
compressors Were carried out at the Naval Postgraduate School during 1983 and 19584. 
It was concluded at that time that the code was able to predict the general behaviour 
of the flowfield with certain discrepancies [Ref. 4]. However, the experimental data in 
each application were not complete enough to enable the evalution of each aspect of the 
code independently, such as the axisymmetry assumption in modelling the physical flow 
and the use of two dimensional cascade data for correlating the losses across the blade 
row. Furthermore, the spanwise static pressure distribution, which is the most sensitive 
indicator of accuracy in applying throughflow codes, was not available from the 
experiments. 

In 1985, under the sponsorship of the Naval Air System Command, completely 
detailed flowfield survevs obtained in a five feet diameter two stage axial flow 
compressor at United Technologies Research Center (UTRC) were used to assess the 
accuracy and utility of throughflow codes for predicting compressor flow fields and 
compressor performance. The relevant parameters of the flowfield required for a 
comparison with code predictions, such as the spanwise variation of the total pressure 
loss coefficient, flow turning angle and tangential aerodynamic blockage, were obtained 
from the surveys. Dring and Joslyn [Ref. 5: pp. 56-64] used the reduced experimental 
data in a throughflow finite element code developed by Habashi and Youngson. They 


concluded that an accurate description of the tangential aerodynamic blockage, resulting 


from the non-axisymmetry of the phvsical flow and the presence of wakes, was essential 
if the flow field was to be predicted accurately. The detailed set of measurements was 
subsequently made available as a benchmark data set for the assessment of other 
throughflow codes in reference 5. 

With the complete survey data available as inputs, the Q3DFLO’81 code can be run 
on the benchmark compressor without accessing the two dimensional cascade loss and 
deviation correlations and the endwall boundary laver subroutines. This enables the 
effects of blade wakes and accuracy of the axisymmetric flow assumption to be assessed 
by running the code with and without accounting for the tangential aerodynamic 
blockage. In the process, the accuracy of the code’s numerical method can be examined. 

The benchmark axial flow compressor had extensive hub corner stall in the second 
stage stator. Although cascade data are not appropriate for correlating the losses and 
exit deviation angle across the blade row of compressor with extensive region of stall, it 
is interesting to see how different, quantitively, the predictions will be using the 
correlation subroutine in Q3DFLO’81, which 1s based on two dimensional cascade data, 
compared to the use of actual survey data. Thus interest arises because Dring has 


observed the occurrence of hub corner stall in all the configurations tested at UTRC. 


B. BLOCKAGE FACTORS IN THE Q3DFLO’8!1 FORMULATION 

There are 2 tvpes of blockage involved in throughflow modelling; namely, endwall 
blockage and tangential (two dimensional) blockage [Refs. 2,6]. Endwall blockage 
accounts for the reduction in throughflow area due to the displacement thickness of the 
boundary laver along the shroud and hub surfaces. Q3DFLO’81 accounts for the effect 


of endwall blockage by applying the endwall boundary conditions 
Y= —61(0V nN nue 
and 


m x 
= 52 oF 0» (p an era 
at the hub and case walls respectively. Here, m is the massflow rate and oy and 63 are 
the boundary laver displacement thicknesses at the hub and shroud respectively. 
The tangential blockage, (1 - b), accounts for the reduction in the effective annulus 
area in the circumferential direction. The tangential blockage varies as a function of 


radius. Dring [Ref. 5: pp. 11] states that “the blockage (1 - b) reflects the fraction of the 


10 


at 


circumference not available to the inviscid throughflow.” Tangential blockage can be 
mechanical or aerodynamic. Mechanical blockage occurs within the blade row due to 
the airfoil tangential thickness. Aerodynamic blockage occurs both within the blade row 
as well as in the annulus region between blade rows. Whithin the blade row, the 
aerodynamic blockage results from the displacement thicknesses of the boundarv lavers 
along the blade surfaces and the departure from axisvmmetry. In the annulus region, 
the aerodynamic blockage is due to blade wakes and non-axisvmmetries such as corner 
Stall. 

Q3DFLO’8! accounts for the tangential blockage bv defining the axial (1¥’,) and 


radial (IV,) velocity components in its governing equation as 
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However, Q3DFLO’8! only computes the tangential blockage (1 - b) due to the 
mechanical blockage within the blade rows, and does not include aerodvnamic blockage. 


There 1s no provision for introducing the aerodynamic blockage as input data. 


C. ASSESSMENT METHODOLOGY 

Measured values of deviation angle and loss coefficients can be used as inputs to the 
code. Consequently, the applicabilitv of the correlation expressions for the losses and 
deviation angle that are contained in Q3DFLO’S1, to the prediction of the performance 
of a compressor with a large stall region, can be assessed. First, the code can be run 
inputting measured quantities from the benchmark data set, and then run again using 
the correlation subroutines. Both runs will use the same geometrical data and will have 
the same inlet conditions. 

An assessment of the effect of including tangential aerodynamic blockage cannot 
be done in a straight forward manner since Q3DFLO’81 does not have provisions for 
introducing tangential aerodynamic blockage as input data. What can be done 1s to use 
an average value for the tangential aerodynamic blockage (given by reference 5 on page 
21 as a mass-average value) and equate it to the equivalent endwall blockage, which the 


code can accept as an input. 


This approach, although not physically correct, does give an indication of the effect 
of the tangential aerodynamic blockage on the throughflow prediction. In a broad 
sense, both types of blockages are similar in their effect. They both increase the 
meriodional velocity component in the throughflow. From the results presented by 
Dring and Joslyn [Ref. 5: pp. 49,99,102], one can conclude that, in the absence of a 
detailed description of the radial variation in the tangential aerodynamic blockage, a 
uniform tangential aerodynamic blockage (an average value) produces much better 
results than bv not including the tangential aerodynamic blockage at all. 

Three cases were run on the benchmark compressor using Q3DFLO’81. Cases 1 and 
2 were run using the measured values of loss and deviation angle. However, Case 1 was 
run with no tangential aerodynamic blockage and Case 2 with an equivalent tangential 
aerodynamic blockage added to the case-wall. The results of Cases 1 and 2 are 
compared with the benchmark compressor experimental results to assess the effect of 
including tangential aerodynamic blockage on the prediction. Case 3 was run using the 
Q3DFLO’S! correlation subroutines and with the same equivalent tangential 
aerodvnamic blockage added to the case-wall as for Case 2. The Case 3 prediction was 
compared to those of Cases 1 and 2. In all the three cases, the endwall boundary layer 


computation in the code was not used. 


D. INPUT DATA 

The annulus of the benchmark compressor was cvlindrical. The radius of the hub 
and tip were 24 ins and 30 ins respectively. The airfoils of the rotor and stator were 
NACA 65 series with circular arc camber lines. Figure 5 and Figure 6 show a schematic 
and the computational mesh for the benchmark compressor respectively. Only the 
second stage rotor and stator were modelled. This was because the detailed 
measurements Were provided for the second stage and furthermore, Q3DFLO’81 1s 
limited by the grid capacity and could not model the whole compressor and maintain 
adequate resolution near the walls. A total of 15 stations was used. Station 1, the inlet 
station to the second stage, is located mid-way between the first stage stator trailing edge 
and the second stage rotor leading edge. Stations 2 and 6 correspond to the second stage 
rotor leading and trailing edges respectively, and Stations 8 and 12 to the second stage 
stator leading and trailing edges respectively. The code predictions and benchmark 
measurements were compared at Stations 7 and 13, which correspond to Stations 4 and 


5 in reference 5. 


A total of 11 axial grid lines (which is the maximum permitted by Q3DFLO’81) was 
used. The streamlines were also clustered towards the shroud by using a repartition 
factor of 0.1. This emphasis on the tip region was necessitated by the abrupt variation 
in the benchmark loss coefhcient near the casewall. Figure 7 shows the distributions 
of the velocities and total pressure as a function of radius at Station 1, which were input 
for all three cases. Figure 8 and Figure 9 show the distributions of the benchmark loss 
coefficient and exit flow angle as a function of radius for the second stage rotor and 
stator, which were input for Cases | and 2. 

Several points are noted. The benchmark data are given in terms of 
non-dimensional quantities. For use in Q3DFLO’81, the data were converted to the 
dimensional form in metric units. A Standard Day of 15° C and 1.01325 bar was chosen 
as the inlet condition of the compressor to calculate the mass flow rate and the 
velocities. The loss coefficients given in the benchmark data are specified with respect 
to the mid-span wheel speed. For Q3DFLO’8], the loss coefficients were re-calculated 
with respect to the benchmark rotor inlet relative velocity and the stator inlet absolute 
Velocity at their specified radial positions, for the rotor and stator respectively. In 
Q3DFLO’81. the aerodynamic input data are limited to 10 radial points for each 
variable. Since I$ to 20 points were specified in the benchmark data set for each 
variable, ten points were carefully selected such that thev represented the macro 


variation as closely as possible. 


E. EFFECT OF TANGENTIAL AERODYNAMIC BLOCKAGE ON THE 
THROUGHFLOW PREDICTION 

Figure 10 through Figure 17 compare the Q3DFLO’81 code predictions for Cases 
1 and 2 with the experimental results. The comparisons are made at Stations 7 and 13 


for the following parameters: 
e Total pressure coefficient 
e Relative total pressure coefficient 
e Static pressure coefficient 
¢ Axial velocity component 
e Absolute flow velocity 
e Relative flow velocity 
e Relative flow yaw angle 


e Absolute flow yaw angle 
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Q3DFLO’$8!1 outputs the pressure and flow velocitv in metric units; namely, bar, and 
meter per second. These parameters, shown plotted in Figure 10 to Figure 17, were 
non-dimensionalised using the mid-span wheel speed. 

At Station 7, Figure 10 through Figure 17 show that the predictions of Cases | and 
2 are verv similar. It can be seen that Case 2 agrees slightly better than Case | with the 
benchmark data in both magnitudes and overall trends. 

At Station 13, Figure 12, Figure 13 and Figure 16, however show that there are 
significant differences in the predictions of static pressure, axial velocity and relative flow 
yaw angle between Case 1 and Case 2. The prediction of Case 2 agrees very wellauarna 
the benchmark data, but not Case 1. Case 1 shows appreciably higher static pressure 
and relative flow yaw angle and lower axial velocity. 

The agreement at Station 7 and disagreement at Station 13 between Case | and Case 
2 are, almost certainly, the result of the tangential aerodynamic blockage introduced at 
these stations in Case 2. The tangential aerodynamic blockage 1s 1.4 percent at Station 
7 compared to 4.7 percent at Station 13. Fora constant massflow rate. the absence of 
1.4 percent area blockage at station 7 in Case | causes a slight reduction in the axial 
velocity Which results in onlv a slight increase in static pressure and absolute flow vaw 
angle as can be seen in Figure 12 and Figure 17. However, the effect of 1snoimmemune 
4.7 percent blockage at Station 13 in Case 1 1s clearly significant. It causes a marked 
reduction in the predicted axial velocity, which results in a significant increase in the 


Static pressure and relative flow vaw angle. 


F. CODE PREDICTION USING TWO DIMENSIONAL CASCADE DATA 

Figure 10 through Figure 17 also show the predictions for Case 3, which uses the 
code’s correlations of cascade data, rather than measured data, for computing loss and 
deviation angle. The Q3DFLO’$! endwall boundary layer part of the code was not used 
because of the limited resolution in the benchmark data and difficulties in starting the 
endwall boundary laver calculation in the middle of the compressor. Clearly, the most 
meaningful evaluation of the accuracy of the cascade correlation will come from a 
comparison of measurements and predictions in the core flow, say from 25 percent to 
75 percent span from the hub, 

Within the core flow, the Case 3 predictions of the velocities and flow angles at both 
Stations 7 and 13 generally follow the trends in the benchmark experimental results. 
However, they are not as good as the predictions obtained in Cases | and 2. The 


maximum disagreement in the flow angle is seen to be less than three degrees and in the 
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flow velocity, less than ten percent. The code also computed much higher total and 
static pressures, especially at Station 13. This would be expected perhaps, because of the 


higher losses resulting from the presence of hub corner stall. 


G. OBSERVATIONS 

Comparison of the Q3DFLO’81 predictions for Cases 1 and 2 clearly demonstrate 
the need to include the tangential aerodynamic blockage resulting from the wakes, or 
an equivalent measure of the departure from axisymmetry, in computing the compressor 
throughflow. The good agreement of the predictions of Case 2 with the benchmark 
experimental results, is considered to validate the code’s numerical method. 

Results from Case 3 also show that it is necessary to know in advance whether the 
compressor being modelled has an extensive stalled region of flow. If the compressor 
has an extensive region which Is stalied, the subroutines contained in Q3DFLO’S1 for 
calculating the total pressure losses and deviation angles will not give accurate results. 

However, in general, given the correct inputs, the s1mplified throughflow model 1s 
able to predict correctly. both the levels and the spanwise distributions of the flow at 


each station 1n an axial flow compressor. 
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IV. APPLICATION TO A MULTI-STAGE COMPRESSOR 


A. BACKGROUND 

The three-stage, axial flow compressor test rig at the Turbopropulsion Laboratory 
of the Naval Postgraduate School (NPS) 1s presently configured for an experimental 
program to determine the effect of rotor tip clearance on the efficiency and flow behavior 
in an axial flow machine. Figure 18 shows a radial section of the compressor. It has 
three repeating stages of symmetrical blading. It also has a row of inlet guide vanes to 
provide the required pre-rotation to the first stage rotor. The stator provides the same 
pre-rotation to each following stage. The design, geometry and construction of the 
compressor are described in references 7 and 8. 

Q3DFLO‘’S1 was applied to analyse the machine performance and flow through the 
test rig. The purpose was to see whether Q3DFLO’81 could be used to support the 
experinental research activities on the effect of rotor tip clearance. Furthermore, if 
Q3DFLO’81 was shown to predict the performance of the test rig correctly, then the 
predicted throughflow solution could be used as input for the blade-to-blade analysis, 
which ts also available in Q3DFLO’S81. 

The present chapter deals with the application of the throughflow code contained 
in Q3DFLO’S1 to the NPS compressor test rig. In the process, some practical aspects 
of throughflow modelling using Q3DFLO’81 were uncovered. The configuration of the 
test rig modelled here consisted of the IGV’s and only the second and third stages, since 
this was the configuration being used in the experimental investigation of the effects of 


tip clearance change. 


B. EXPERIMENTAL SET-UP 

Figure 18 also shows the experimental arrangement that was used to measure the 
overall stage pressure rise of the compressor. The overall static pressure rise was 
measured using two static pressure taps located on the shroud 1.6 inches ahead of the 
second stage rotor and two taps at the exit of the third stage stator. Each pair of taps 
were spaced about 90° apart. 

The mass flow rate was calculated from an inlet bellmouth measurement, involving 
four impact pressure probes and four static pressure taps spaced equally apart, and one 


thermocouple sensor. 
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For the purpose of prescribing the inlet conditions for the code, the flow at the inlet 
to the [GV was measured. The total temperature was assumed to be uniform and the 
same as for the bellmouth measurement. The spanwise distribution of the total pressure 
was measured using a nine probe total pressure rake located 0.87 inches downstream of 
the exit plane of the strut which was 1.75 inches upstream of the IGV. The nine probes 
were aligned along a radial line positioned almost mid-way between two adjacent struts. 
The two probes nearest the endwalls were 0.2 inches and 0.25 inches away from the hub 
and shroud respectively. The static pressure was measured at the shroud and the hub 
using two static taps, and at mid-channel using a pitot static pressure probe. The 
Velocity was assumed to be axial and its magnitude was calculated from the dynamic 
pressure distribution. 

All measurements Were recorded using a Hewlett Packard 3053 Data Acquisition 
System. A description of the data acquisition system, the instrumentation and the 
measurement uncertainties are given in reference 9. The uncertainty of pressure rise 
coefficient and flow coefficient was established from repeated measurements to be 0.5 


Pereent and 0.6 percent respectively. 


C. EXPERIMENTAL MEASUREMENTS 
Experimental measurements were obtained at a fixed compressor speed for four 
mass flow rates. For performance mapping, the mass flow rate is expressed in terms of 


the dimensionless flow coefficient (¢) defined as 


1 
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where 7m is the massflow rate determined from the bellmouth measurement, p,,, and 4,,, 
are the density and annulus area respectively at the IGV inlet, and U, is the mid-span 
peripheral speed of the rotor. 

The flow coefficients so obtained were 0.61, 0.67, 0.70 and 0.75. Table | tabulates 
the pertinent experimental measurements and calculated parameters. The static pressure 
rise (AP,) across the rotors and stators is expressed by the dimensionless quantity, IT,, 
which is defined as 

AP, 
= 


] 2 
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The variation of I], with @ is the compressor ‘characteristic’. In the following section, 


the measured characteristic is compared with the Q3DFLO’81 code prediction. 


Table 1. EXPERIMENTAL FLOW AND PRESSURE COEFFICIENTS 
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The measured spanwise distributions of total and static pressure and the velocity 
calculated from them. are shown in Figure 19. The static pressure 1s observed to be 
skewed linearly from hub to shroud based on the shroud, mid-span and hub 


measurements. The (axial) velocity was calculated from the measurements using 


Vy = lucy iov ee 


a Pigy 


where the density (,,) in pounds mass per cubic foot at the IGV inlet was obtained 


using 


— _ Psigv 518.69 
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x 0.07647 
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D. MODELLING OF THE COMPRESSOR 

The NPS 3-stage axial compressor has a cylindrical annulus. The radius of the hub 
and tip are 10.8 inches and 18 inches respectively. The airfoils of the rotor and stator 
are both of the C series type with a modified thickness distribution. The compressor 
was modelled from the measurement plane upstream of the JGV to the station 
downstream of the third stage stator. Figure 20 shows the computational mesh, which 
had a total of 25 stations and 10 axial grid lines. Grid lines were clustered towards the 
shroud by specifying a repartition factor of 0.1. As the total number of nodes then 
exceeded the maximum of 600 nodes allowed in the Q3DFLO’81 code, the analvsis was 
done in two parts. The first part analysed the flow from Station 1, ahead of the IGV, 
to a station immediately downstream of the third stage rotor. The prediction at this 
station was then used as input to a second calculation from the inlet of the third stage 
Stator to the measurement station downstream of the stator exit. 

The input data required by the code for the inlet condition at Station 1 were the 
mass flow rate, the spanwise distributions of the velocity components. the total 


temperature and the total pressure. 


FE. EFFECT OF THE STRUT WAKES 

The boundary layer displacement thickness at the shroud was measured to be about 
Oeleeiniches in a separate study by Tarigan (Ref. 10]. For the code prediction, a nominal 
endwall blockage of 2.6 percent was specified at Station 1, which 1s equivalent (in area) 
to a 0.15 inches displacement thickness at the shroud. The higher blockage factor was 
to account for the thicker boundarv laver at the junction of the strut and the endwalls 
and the negligibly thin boundarv laver at the hub. 

A disagreement with the mass flow measurement obtained from the inlet bellmouth 
was noted when the spanwise velocity distribution at the IGV inlet was integrated to 
determine the mass flow rate. The data in Table 2 show that the massflow rates 
determined from the measurements at Station 1 were about 3.5 percent higher than the 


bellmouth measurement for all four flow coefficients. 
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Table 2. BELLMOUTH VS INLET RAKE MEASUREMENTS OF FLOW 
RATE 
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Examination revealed that the measurement plane at Station | was about 21 percent 
of chord behind the strut trailing edge. There were a total of six struts equally spaced. 
Figure 2! shows the profile of the strut which had a chord of 4.12 inches and maximum 
thickness of 0.813 inches. It was known that strut wakes were present at the 
measurement plane, although they did not intercept the rake sensors. The thickness of 
the wake had not been measured since the mechanical arrangement of the casing did not 
permit a circumferential survey of the flow. However, two static pressure taps on the 
shroud, located between the IGV outlet and the 2nd stage rotor inlet, did indicate an 
increase in the static pressure, as the flow proceeded downstream. This could be 
explained by the diffusion of the strut wakes in the annulus. The presence of the strut 
wakes also explains the apparent increase in mass flow rate mdicated by the rake at 
Station |. The strut wakes, were effectively a two dimensional blockage. causing a 
reduction in the circumferential flow area and an increase in the local flow velocity from 


hub to shroud at the IGV measurement plane. 


F. ADAPTATION OF THE INLET FLOW CONDITION 

The displacement thickness of the strut wakes could be expected to be fairly uniform 
in the core of the flow but not at the corners between the strut and the endwalls. For 
numerical modelling. it was reasonable to assume that it was uniform from hub to 
shroud. The endwall boundary layer displacement thickness, the mass flow rate 
determined from the bellmouth measurement and the spanwise axial velocity measured 
at Station | were sufficient to calculate a unique displacement thickness of the strut 
wake which gave the measured mass flow rate. 

The wake displacement thickness can be viewed as tangential blockage. 
Unfortunately, the Q3DFLO’81 code was unable to accept a tangential blockage factor 
as input data. Consequently, it was assumed that the strut wakes were completely 


diffused at constant radi before entering the IGV, and this diffusion process was 
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calculated externally from the Q3DFLO’81 code. Appendix (B) describes the diffusion 
calculation for the strut wakes. The spanwise distributions of the diffused (or ‘mixed 
out’) axial velocity and total pressure were then used as boundary conditions for the 
code, at Station |. 

One could expect that, by assuming the wakes had diffused before entering the IGV, 
instead of letting them be swallowed by the IGV’s and then diffused gradually to have 
negligible thickness at Station 8, the prediction of the axial velocity component would 
differ from the actual value immediately downstream of the IGV. But at Station 8, and 
beyond where it was expected that the strut wakes would have diffused completely, the 
code prediction would be valid. A comparison of the code predictions with the 
experimental measurements for the static pressure at the shroud at Station 8 would give 


a good indication of the validity of the above reasoning. 


G. RESULTS FOR THE ADAPTED INLET FLOW 

The endwall blockage and wake thickness of each strut at Station | required for the 
code calculation to match the measured static pressure at Station 8, are shown in 
ole: 3. 


Table 3. DISPLACEMENT AND WAKE THICKNESSES AT STATION 1 
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It was found that a wake displacement thickness of 0.528 inches and an endwall 
blockage factor of 0.976 produced a match with the measured static pressure at Station 
8 fairly well. The ratio of the wake displacement thickness to the strut maximum 
thickness was 0.59. 

Figure 22 shows the spanwise distributions of total and static pressures and axial 
velocity after the diffusion. It can be seen that the total pressure was hardlv affected. 
AS expected, there was an appreciable increase in the static pressure of about 0.5 inches 


O@iwater dile tela 3.5 percent reduction in velocity. 
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H. COMPARISON OF THE OVERALL PERFORMANCE 
Q3DFLO’S1 was used to predict the overall performance at each of the four flow 

coefficients both with and without strut wake diffusion at Station 1. The computed 
variation of FI, vs @ is shown in Figure 23 denoted as follows: 

e Case 1]: Experimental Measurement. 

e Case 2: Q3DFLO’81 Prediction with Correction for Strut Wake Blockage. 

¢ Case 3: Q3DFLO’81 Prediction without Correction for Strut Wake Blockage. 

e Case 4: Meanline Analysis {Ref. 9] 


It should be noted that the flow coefficient of Case 3 was based on the actual mass flow 
rate rather than the apparently higher mass flow rate measured at the IGV inlet. 

Comparing Cases 2 and 3, Case 2 shows significantly higher pressure rise coefficients 
than does Case 3. This is probably for two reasons. First, the reduced axial velocity 
Which results from wake diffusion causes an increase in the blade loading. The IGV exit 
flow angle calculated by the code 1s independent of the value of the flow coefficient. The 
code predicted essentially the same absolute flow angles at all stations from the exit of 
the IGV at Station 4 to the inlet of the first rotor at Station 9, both for Cases 2 and 3. 
With the same absolute flow angle, and lower axial velocity, Case 2 involves a 
correspondingly lower absolute tangential velocity upstream of the rotor. This would 
cause an increase in the relative tangential velocity component, which results m an 
increase in the blade loading of the first rotor. Second, the total pressure loss coefficient 
calculated by the Q3DFLO’S1 subroutine for Case 2 was significantly lower than for 
Case 3 for all blade rows. This suggested that Case 2 involved more favorable incidence 
angles than Case 3. 

Comparing Case 2 with Cases 1 and 4, Case 2 showed better agreement with Case 
1 (measurements) at the flow coefficient of 0.612 but departed by an increasing margin 
for flow coefficient above 0.67. Since the design flow coefficient was about 0.64, the 
prediction was better near the design condition than at off-design conditions, which 


might be expected. 


I EFFECT OF THE IGV EXIT FLOW UNDERTURNING 
It was suspected that the difference in the pressure rise coefficient between the 
experimental measurements and the predictions at all conditions might be the result of 


underturning of the IGV exit flow, which had been observed in the compressor. 


22 


Waddell [Ref. 11! had measured the IGV exit flow angle and had established that 
it was independent of the flow coefficient. He showed that the turning angle was smaller 
than the design intent over a large part of the span. The Q3DFLO’81 prediction of the 
IGV exit flow angle for Case 2 1s compared with Waddell’s measurement in Figure 24. 
It can be seen that the flow was underturned over 70 percent of the span and the 
maximum underturning was about five degrees. 

Since the metal trailing edge angle of the IGV could be adjusted readily in 
Q3DFLO’81, attempts were made to simulate the measured IGV exit flow angle in Case 
2 by adjusting the IGV metal angle. This case is denoted as 

e Case 5: Case 2 with Simulated Measured IGV Exit Flow Angle 


Cases | and 2 are shown replotted in Figure 25 with the prediction for Case 5. Case 
5 shows a fairly constant improvement over the measured range of flow coefficients. It 
reduces the difference between prediction and measurement of the pressure rise 
coefficient by about 40 percent at the on-design flow coefficient but has progressively 


less effect for higher flow coefficients. 


J. OBSERVATIONS 

From the above studies, it is clear that it is desirable to have provision for 
introducing tangential aerodynamic blockage as input data. The presence of the struts 
ahead of the IGYV is typical for an axial flow compressor. Its effect on the flow velocity 
at the inlet of the [GV must be considered, and this can be accounted for by using the 
tangential aerodynamic blockage factor. 

The prediction of the performance near the design condition 1s much better than at 
off-design conditions. An accurate description of the flow conditions at the inlet of the 
compressor is essential if the code is to give good predictions at following stations. 

The prediction of the IGV exit flow angles bv Q3DFLO’81 differed from the 
measured values. As a result, the predicted flow downstream of the IGV exit was not 
representative of the physical flow. It is noted that the code does not contain a separate 
correlation for deviation angle from IGV’s, but treats the IGV’s as a stator (compressor) 


blade row. 
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V. THE INVISCID BLADE-TO-BLADE SOLUTION 


A. BACKGROUND 

Q3DFLO’81 has a separate inviscid code to solve for the velocity and pressure fields 
on the axisymmetric streamsurface between two adjacent blades. The blade-to-blade 
code requires the inlet and outlet flow conditions to be specified. The code can be run 
independently by entering the inlet and outlet conditions manually. Alternatively, 
Q3DFLO’$1 can optionally couple the throughflow code with the blade-to-blade code. 
In this case, Q3DFLO’81 would first compute the results for the throughflow. It then 
transfers the necessary data internally to the blade-to-blade code, for the specified blade 
row, and continues with the blade-to-blade computation. 

The blade-to-blade code was apphed to Cases 2 and 5 (of the preceding chapter) for 
flow coeflicients of 0.61 and 0.67 and the inviscid blade-to-blade solution for the second 


rotor at the tip section was generated. 


B. INPUT DATA AND SOLUTION 

The blade profile and the computational mesh for the rotor up calculation are 
shown in Figure 26 and Figure 27 respectively. A total of 25 stations was used. 
Stations 6 and 20 are the inlet and outlet stations of the rotor blade passage. There are 
two options for specifving the outlet flow condition. Either the Kutta condition 1s 
imposed at the trailing edge or the exit relative flow angle must be specified. The 
solution presented herein is based on the use of the Kutta condition. 

As an example of the blade-to-blade computation, the results are presented here for 
Case 5 at a flow coefficient of 0.61. The computed iso-pressure lines on the 
axisymmetric stream surface at the tip and the pressure coefficients on the suction and 


pressure sides of the blade, are shown in Figure 28 and Figure 29 respectively. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


The finite element code Q3DFLO’S1 was applied successfully to two different axial 
compressor geometries. First, application to a ‘benchmark’ compressor, for which very 
complete flow survey data were available, showed the importance of accounting for 
peripherally non-uniform flow (blockage) effects. Second, application to the NPS axial 
research compressor demonstrated clearly the difficulties that are inherent in comparing 
code predictions with experimental data when the experimental information is not fully 
comprehensive, and when the measured flow is not valid as an inlet boundary condition 
for tne code. 

In spite of the complex nature of the flow within the NPS compressor test ng, when 
a valid inlet boundary condition was derived from the measured data, the throughflow 
code of Q3DFLO’S1 predicted the pressure rise fairly well near the design operating 
condition. An analvsis of the measurements, in order to derive a proper boundary 
condition for the code, was required in order to achieve agreement. At off-design 
conditions, the code was not as successful. For flow near to the design condition, the 
inviscid blade-to-blade code was used successfully to generate the 1so-pressure lines and 
the blade surface pressure distributions for the tip section of the second rotor. 

From the above experience, the following were found to constrain the application 
pe OSDFLO'S1: 


e The code does not accept a tangential aerodvnamiuc blockage factor as an input for 
the throughflow prediction. This did not allow the strut wakes at the IGV inlet 
measureinent plane to be introduced to the throughflow in the manner that they 
occurred physically. 


e The code allots the endwall blockage factor at the inlet station equally between the 
hub and the case walls. This was not representative of the flow in the compressor 
test ng. The blockage at the inlet was largely due to the boundarv laver 
displacement thickness at the case wall. The hub had a negligibly thin boundary 
knoe 


e The code only permits the meriodional mesh lines to be clustered either towards the 
case wall or the hub. In applving Q3DFLO‘’81 to the benchmark data set and the 
compressor test rig. it was found that it would be desirable if the meriodional mesh 
lines could be clustered towards both the case wall and the hub. This is because 
abrupt changes in the flow conditions tend to occur at both the case wall and the 
hub. 


e The code limits the capacity of the mesh to 600 nodes and the input values for each 
variable such as the total pressure at the inlet station to 10 points. A larger 
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capacity is required to model a multi-stage geometry using grid lines clustered near 
the walls. 


e The use of two dimensional cascade data for calculating losses across the rotor and 
stator blade rows was shown not to be accurate for a blade row with corner stall. 
It is therefore advisable to establish whether stall regions are present in any new 
compressor to which the code 1s applied. 


¢ The code did not predict the exit flow angle of the inlet guide vane correctly. 
Although this was overcome by adjusting the geometry of the inlet guide vane in 
the code, the correction required that the actual flow angle be known. This 
information would not generally be available in a normal application of the code. 

The use of the code in the tip-clearance investigation on the NPS compressor is seen 
to be limited by the above constraints, but also by constraints on the experiment itself. 
For example, the spanwise distribution of the strut wake profile at the inlet measurement 
plane needs to be measured accurately. Although the assumptions of uniform wake 
thickness and wake diffusion on the constant radi axisymmetric stream surface worked 
fairly well, an accurate representation of the inlet velocity near the case wall 1s required 
to obtain an accurate prediction of the pressure rise, and a more precise inviscid solution 
at the rotor tp. 

Unfortunately, the compressor is not in general, equipped for peripheral flow 
surveys, and the required modification to the heavy case wall would be difficult and 
expensive to make. Jt 1s clear however, that peripheral survey data and the means to 
input such data into the code, are essential. 

Finally, it is noted that the present assessment was made using an early edition of 


the code. Later versions of the code might be less constrained. 
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Figure 1. Stream Surfaces Through an Annular Blade Row 
Source: Varva M. H., Aderothermodynamics and Flow in Turbomachines, Wiley, New 
York, May 1974. 
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Iso-pressure Lines from the Blade-to-Blade Solution 
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Figure 29. Pressure Coefficieuts from the Blade-to-Blade Solution 
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APPENDIX A. DERIVATION OF THE RELATIVE THROUGHFLOW 


EQUATION FOR Q3DFLO‘81 


The steady, relative equation of motion (Eq. 18), is 


—_ 


W x (Vx V) = TVpS —Val + + Fy 


{18} 


The two force terms F, and F, are eliminated from Eq. (18) ) by projecting Eq. (18) in the 


direction of F, x W. Since F, iS Opposite in direction to W 
(F, x W)+F,=0 
and, clearly 
(F, x 1). F, =0 
Hence, the balance 0: Eq. (13) become: 


—(F, x W)- Wx (Vx VF) =(F,x 1) «(TV 2S — Val) 


Next, Eq. (22) 1s expressed in cylindrical coordinates. [his is doneterm by tem 


Step 1: Expanding F, sti gives 
Fy x W = i,(Fyg ll, — Fogo) + ig Foo, — Fop!¥_) + ig( Fo Vo — Fog lV) 
Step 2: 7V,S —V,/ on the RHS of Equ(22 ibecomes 


Cl = oS Cl = A 
Nia) as Ora 70 — 39 ) tills — 











eS 
TV pS — Vpl = i,( Ta — 


on) 
~ 


{20} 


{21} 


{22} 


{23} 


{24} 


Since the model assumes axisymmetric flow by averaging the flow properties across the 


blade passage 


1 o..) | 
: CO 





{25} 


and Eq. (24) is reduced to 


a bea a 


Sam! ey oS. OT 
ns — Val = i,(T— >) rae ae) {26} 


Sep. Using Eq: (23) and Eq. (26), the RHS of Eq. (22) becomes 


= ae i 
(Fy x W)°(TVpS — Vel) = (Fig WV, — Fz Wo) 


OQ) 






































Or ) 
CS 27) 
f C C 
Step 4: Since 
= =) oY, ari). = a¥, av, =) av) ar 
iene Sense Teen ace es a ee 8 
meme Eq. (25) 
ee celia OM or = 1 Sr’) x 
Wea —— l, ; a: @ az aa Gr J+, fe Gr {29} 
therefore 
— i -—— Or Crh g) OV, OVs 
i (van) — =) ae 
=f -—*)) 
— Wy, €(rV6) WV. O(rVe) 
eee 30 
7A ae As : Ar ) {30} 
—- ov, OV W, e(rV,) 
1 ile ~ Sa a 


Since there 1s no difference in magnitude and direction between V, and I]’, and J’, and 
H”,, Eq. (30) can also be expressed as 


s ee ar AW. BI 
Tea ey SOE 











ra 
r Cr ales Gz ~—soOr )) 
Smee Ft) W, c(rVo) 
+(e a ) {31} 
— ay CW Vo a(rVo) 
oS. 


37) 


Step 3; Using Eq: (23) and Eqe(31); theweias On eam seccomnes 


ete Po m W, Ors) OW, aw, 
(Fi, x) «WY! (VY) (Egg Ys a i) ee 


Arn) WY, erp) 














+ (FogI¥, — Fo W Met ES {32} 
CH a cle Wy O(rVy) 
+ (Fo Wo — Fog VVC a — = ) — a a) 


Re-arranging, 














ee — al, ow, os 2 
G(r 
uae 2 (Ee, — Fea, — Fee 
a 
-— ec (F,,1¥> — F100, + F358, — FP 
33 
aw. al, a eS) 
= (= oe) = ate 





TOU) ae 
$F (WFeg Vo + FeeW,) — Fy Vg + 12) 





Or 
Le et a) , ; 
- (ae (Fa( IV, + 1VG) — Wi Fog, + Fog H¥'9)) 


Substituting }]? = 117? + 1V3+ 1? into Eq. (33) yields 








ale aH, 
(F, xa» 1 x (Ve a (Eu — WE gi + Ftv, + Pee 
a a) 
+ = —— (WV Fg Wg + Fp WV, + FoI, — Fy’) 


ry 














| a @) 
a = (F,,W? amet Hy Adeug bee + f,,|V, + Fig lV 9)) 
34 
(oe Me thE 7 
= ( Ap — ae ) be = al oe r) 
C(r} => = 

+r e 2 (HF, » HW) — Fyz1V*) 

joe FB (F,,V° — W (F, . 1V)) 
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Since F, + 1’ = 0 in Eq. (34), the LHS of Eq. (22) becomes 


oot 


(F,x W)- Wx (Vx V)= 





pope St) _ ee OrVo) 1 For OlrVo) (35} 
ji ° 02 f Fog Or “ ng OZ 


eieeeo. From Eq. (27) and Eq. (35), Eq. (22) becomes 























ge oe rN) Fe 1 KS) 
or Oz IO auaae 6 ee lz 
{36} 
=—> ((W,- Po; WoT aS eReC! )+( oor yy W \(T eS mea, 
(a nO a Gz 


Step 7: The body force components from Eq. (36) can be eliminated by expressing them 


in terms of the blade surface angles. Figure 30 shows that, 


oe 
—* = —tan p (ou 
Fog 
and 
jee | 
a — Coney, 38 
ee {38} 














Cr Oz 
A és ee cS 
W’, + W, tan B)\( —— - TS ; —W - 
esa eae foe) Uh, tany — I,)(-—-— T——) £39} 
W? 
| o(r V5) C(r Fo) 





caters CONG eeepc UT) 17 
Cr Cz 


Step 8: Lastlv, the streamfunction (W) defined through the equations 


“ 


ey zs = p2arblV, {40} 
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and 


~ 


ow r . 
= ——7/) / 4 
2n Ap p2rrblyV, {41} 





where b is the tangential blockage ratio, is introduced. The LHS of Eq. (39) can be 


expressed in terms of w by differentiating Eq. (40) and Eq. (41), which yields 








cW,  @,_1 oY 
= — (-—— ——— 4 
Cr Gr * prb ¢ér ) a 
and 
OV _ A 
a ae (43) 


Finally, Eq. (19), which is the radial equilibrium equation for the rotor flow, 1s obtained 


by substituting Eq. (42) and Eq. (43) into Eg. (39), yielding 


a ) ew a 1 ¢y 



































Oey ( prb Cr ie GZ ( prb €z )= 
“ed oS ol oS 
Wi + 1, ti — , — W = 
(H+ TV, tan B)( ay Poa, + OF, tan n al A. rx) (44) 
12 
1, C(r V9) e(r Vg) 
cael ale Ay + tan 7 As 
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APPENDIX B. CALCULATION OF THE MIXED-OUT CONDITION FOR 


THE STRUT WAKES 


The process of the diffusion and mixing of the strut wake on an axisymmetric stream 


surface at a constant radial position can be modelled as constant area two-dimensional 


mixing of two streams initially separated by the wake thickness (6%) as shown in 


Figure 31. Since the flow Mach number 1s less than 0.2, the velocity and total pressure 


after mixing can be determined by assuming the flow to be incompressible. 


conservation of mass, 


ox 


PS OO 


so that 





y= (1 es. Vy 
Using the definition of blockage factor (k) 


k=(1-) 





Eq. (46) becomes 
ct 
By conservation of momentum, 
(P, — Py)s + pls— 8) -— ps, =0 
Using Eq. (46), Eq. (49) becomes 
(P, — P,)s + pls — 8°)V2 — psi - Srv? =0 


Re-arranging, 





By 


{45} 


{46} 


{47} 


{48} 


{49} 


{90} 


{91} 


Using Eq. (47), Eq. (S51) becomes 
(P, — P,) = pVi(k — (1 — &)) 

From the definition of total pressure 

1 gies Se > pvt 
and 

ee ie > pV3 
Using Eq. (48), Eq. (54) becomes 

PS + pVek? 
Seetracting Eq. (53) from Eq. (55) 
(Pa — Py) =(P2— Py Fe PVR? = W 
Using Eq. (52), Eq. (56) becomes 
(Pa Pu) = Vik = (1-4) +> obi? = 1) 
Re-arranging, 
Pa = Py — > Vl — 2k +k’) 


Or 


l 
Pa = Py — > eV(1 — ky 


{52} 


{54} 


{97} 


{58} 


{99} 


Hence, the velocity and total pressure after diffusion and mixing of the strut wakes can 


be calculated by using Eq. (46) and Eq. (59) respectively. 
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